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We analyze the magnetization, magnetic torque and susceptibility data of La2- :E Sr a ;Cu 16 ' 18 04 
and YBa 2 63 ' 65 Cu 3 7 -a near T c in terms of the universal 3D-XY scaling relations. It is shown that 
the isotope effect on T c mirrors that on the anisotropy 7. Invoking the generic behavior of 7 the 
doping dependence of the isotope effects on the critical properties, including T c , correlation lengths 
and magnetic penetration depths are traced back to a change of the mobile carrier concentration. 

A lot of measurements of the oxygen and copper isotope effect on the magnetization, susceptibility, magnetic 
torque etc. have been performed on a variety of cuprate superconductors [1-9]. To obtain estimates for the shift of 
the transition temperature T c upon isotope exchange, the respective mean-field models [10,11] have been invoked. For 
the magnetization m they imply that the mean-field transition temperature T c q> can be estimated by extrapolating 
m (T) linearly. To illustrate this traditional treatment we displayed in Fig.l the data of Batlogg et al. [1] for a nearly 
optimally doped Lai.85Sro.i5Cu 16,18 04 powder sample in terms of 16,18 m(T) vs. T. The straight lines are parallel 
linear extrapolations, yielding the mean-field estimates 16 T c o ss 34.56 K, 18 T c o as 34.12 K and AT c q/T c q w —0.012, 
where AY — ( 18 Y — 16 Y) / 1S Y. However, the reliability of these estimates is doubtful because the dominant role of 
critical fluctuations near T c [12-21] and the isotope effects on the penetration depth [22-25] etc. are neglected, the 
anisotropy is not taken into account, and the choice of a linear and parallel portion appears arbitrary. 




FIG. 1. 
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T for a Lai.85Sro.i5Cu 16,18 04 powder sample at 6.5 Oe taken from Batlogg et al. [1]. The 



parallel straight lines are linear extrapolations, yielding the mean- field estimates iB r c0 ~ 34.56 K, 1 
ATco/Tco w —0.012. The black circles are the 16 m (T) data rescaled according to Eq.(4) with a ~ 0.986. 
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Noting that the isotope effects in cuprate superconductors pose a fundamental challenge in its understanding, 
the shortcomings of the traditional interpretation of magnetization and related data call for a treatment that goes 
beyond mean- field and takes the anisotropy into account. In this work we concentrate on the critical regime of 
anisotropic extreme type II superconductors where 3D-XY fluctuations dominate [12-21]. Invoking the universal 
scaling functions for magnetization and magnetic torque we analyze the magnetization data of Lai.85Sro.i5Cu 16 ' 18 04 
[1], magnetic torque data of Lai.92Sro.o8Cu 16,18 04 [6] and the susceptibility data of YBa2 63,65 Cu307-«5 [4] near 
T c . An essential additional relation emerges from the observation that data taken at fixed magnetic field and on 
samples with 16 or ls O, as well as with 63 Cu or 65 Cu, collapse near criticality within experimental error, when 
the temperature is rescaled appropriately. Although this property provides retrospectively partial support for the 
traditional approach [1-9], 3D-XY scaling uncovers the essential role of the anisotropy 7. Indeed, the change of T c 
is found to mirror the shift of the anisotropy 7. As a consequence, the generic shift of the temperature dependent 
magnetization, susceptibility and magnetic torque upon isotope exchange at fixed magnetic field does not provide 
estimates for the change of the transition temperature only, as hitherto assumed [1-9]. Together with the generic 
behavior of the anisotropy [18,26-31], the doping dependence of the isotope effects are then traced back to the change 
Ax u of the underdoped limit x u . It implies a shift of the phase diagram in the temperature-doping plane towards a 



1 



slightly higher dopant concentration x, along with a reduction of the mobile charge carrier concentration 5 = x — x u . 
This contribution leads to a negative shift of T c . We identify a positive shift as well. It stems from the change A70 of 
the critical amplitude 70 at the quantum superconductor to insulator transition. The magnitude and proportion of 
these contributions is controlled by Ax u and A70. Their values are material dependent. In any case, they control the 
isotope effects and remain to be understood microscopically. However, the emerging essential role of the anisotropy 
represents a serious problem for two dimensional models as candidates to explain superconductivity in the cuprates, 
and serves as a constraint on future work towards a complete understanding. 

Whenever 3D-XY fluctuations dominate the magnetization m adopts the scaling form [13-21] 

m(T,6,H) k B Q 3 , 1 dG(z) m 
1/2 

where e (S) = (cos 2 (S) + sin 2 (S) /7 2 ) and 7 = £ a f>/£c denotes the anisotropy. Q3 is a universal constant, G (z) a 
universal function of its argument, £ a b,c the correlation lengths in the a6-plane and along the c-axis, H the mag- 
netic field and <j> the flux quantum. Close to the zero field transition temperature T c the correlation lengths 
diverge as S, a b,c = S,abo.cO \t\ " where v ~ 2/3 and t = T/T c — 1. An essential implication is that in the plot 
m(T,S,H)/ {8)tVh) vs. z = (H^e (S) /$„) 1*1"^ = z~^ v = (*„/ (H^e (S))) 1/2v t the data fall 

close to criticality on a single curve. For YB^CusCv-a this scaling property is experimentally well confirmed [14,16]. 
Because the magnetization exists at T c the combination m (T,S, H) / ^ye 3 / 2 (5) Ty/H^j adopts the universal value 
[15,17-20] 

m(T c ) _ k B Q 3 c 3y00 ^ 



j(T c )e^(S)T c VH $3/2 < 

where C3 i00 is a universal constant [15,17-20]. Thus, plotting m (T) / ^7e 3 / 2 (<5) y/l^tj vs. T, the data taken in different 
fields cross at T c . In powder samples and sufficiently large anisotropy (7 >> 1) this relation reduces to 

m{T c ) 7rfcsQ 3 C3,oo (|cos(5)| 3/2 ^ 



3/2 



1 {T C )T C ^H 2*1 

As the oxygen isotope effect on the magnetization at fixed magnetic field is concerned it implies that the relative shifts 
of magnetization m, anisotropy 7 and T c are not independent but related by Am (T c ) /m(T c ) + A7 (T c ) /-){T C ) + 
AT C /AT C = 0. On that condition it appears impossible to extract T c and AT C from the temperature dependence of 
the magnetization taken in one particular magnetic field. However, there is the special case where close to criticality 
the data scales as 

18 m (T) = 16 m (aT) , (4) 

within experimental error. Subsequently it would also justify the traditional method of extracting AT C /T C [1-9]. 
A glance to Fig.l shows that this scale transformation is well confirmed within experimental error, roughly given 
by the size of the symbols of the data points. It yields Am (T c ) jm (T c ) ~ with a = 18 T C / 16 T C ~ 0.986. Thus, 
AT C /T C ~ -A7 (T c ) /7 (T c ) ~ -0.014 for Lai. 8 5Sro.i 5 Cu 16 ' 18 04. Although AT C /T C ~ -0.014 is close to the traditional 
estimate , AT c q/T c q w —0.012, its significance is fundamentally different. Indeed, because Am(T c ) /m(T c ) ~ 0, the 
universal relation Am (T c ) /m (T c ) + A7 (T c ) /j (T c ) + AT C /AT C = reduces to 

AT C = A 7 (Tc) = AUo A^co = 1 

T c l{T c ) S,abO CcO 

Hence, the isotope effect on T c mirrors that on the anisotropy and the critical amplitudes of the correlation lengths. 
In virtue of the universal scaling expression (2), relation (4) is obtained when the anisotropy 7 scales near T c as 
7 (T) = 07 (aT). The rescaled curves should then cross at T c , provided that the experimental uncertainties do not 
mask the isotope induced change of £ a bo, the critical amplitude of the in-plane correlation length. Considering the data 
shown in Fig.l we observe that the two curves coincide in the critical regime within experimental error, roughly given 
by the size of the symbols of the data points. Apparently, the resolution of the crossing point requires considerably 
more accurate data. Otherwise, as in the present case, the coincidence of the rescaled data confirms the consistency 
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with 3D-XY critical behavior, allows to determine the rescaling factor a around T c , and with Eq.(5) to estimate the 
shifts AT C /T C and A7 (T c ) (T c ) rather accurately. 

Before turning to the implications of these results, revealing that the isotope effect on T c mirrors that on the 
anisotropy 7, it is essential to explore the effect of the dopant concentration. Since sufficiently dense data appears to 
be available for underdoped Lai. 9 2Sr .08Cu 16 ' 18 O4 only, we are left with the reversible magnetic torque data of Hofcr 
et al. [6] shown in Fig. 2 in terms of r vs. T. At T c , fixed orientation and magnitude of the applied field r scales 
as t(T c ) = -constT c -f (T c ) H 3 / 2 [15,17-20] and at fixed magnetic field the shifts are related by At (T c ) /t(T c ) + 
A7 (T c ) h ( T c) + AT C /AT C = 0. From Fig.2 it is seen that with 18 r (T) = 16 r (aT) and a = 0.936 near coincidence 
is achieved within experimental accuracy. Since the universal scaling law for the magnetic torque [15,17-19] is 
essentially analogous to Eq.(2) the near coincidence implies At (T c ) /t (T c ) ~ and (5) holds in this case as well, so 
that AT C /T C ~ -A 7 (T c ) / 7 (T c ) ~ -0.07. 




FIG. 2. Reversible magnetic torque t (6) vs. T at H — 0.1 T and 8 — 45° for La1.92Sro.osCu 
[6]. The black circles are 18 m (T) ~ 16 m (aT) with a = 0.936. 
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To check the generic significance of this scenario we also analyzed the susceptibility data for the copper isotope 
effect in YBa 2 63 ' 65 Cu307_5 of Zhao et al. [4]. For all four dopant concentrations, extending from the underdoped 
to the optimally doped regime, we find that 65 % (T) = 63 x (aT) is satisfied within experimental error, as discussed 
below. Since x = m/H this strongly suggests that the scaling relation (5) holds for both, copper and oxygen isotope 
exchange, irrespective of the doping level. 

Having established the consistency with 3D-XY critical behavior, together with the experimental facts that at 
T c and fixed magnetic field Am/m, Ar/r and A\/x vanish within experimental error, the isotope effect on T c 
does not mirror that of the anisotropy only (Eq.(5)), but is also subject to the other universal relations of the 
3D-XY universality class. In particular, T c , £ c o and A a f,o are not independent but related by the universal relation 
[13,15,19,20,23,32] T c = B^ c0 /\ 2 abQl where B is a universal constant and A Q bo the critical amplitude of the in-plane 
penetration depth. This leads for the respective relative shifts upon isotope exchange to the additional relation 
AT C /T C = A£ c0 /£co — 2AA a {, /A a bo- The lesson is, whenever 3D-XY fluctuations dominate, the isotope effects, e.g. on 
T c ,7 and A a bo are not independent. These relations are particularly useful to open a door towards the understanding 
of the common origin of these effects. For example, given a generic relationship between anisotropy 7 and mobile 
carrier concentration S = x — x u at fixed x, where x u , is the underdoped limit, the isotope effects in the cuprates 
would arise from a shift of x u and the associated change of S. As shown in Fig. 3 for I^-xSr^CuC^, the generic 
doping dependence of 7 is well established in a rich variety of cuprates in terms of the empirical relation [20,31,33] 
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where 70 is material dependent constant. Approaching the underdoped limit, where the cuprates correspond to an 
independent stack of sheets with thickness d s [20,31,33], this relation follows from the doping dependence of the critical 
amplitudes of the correlation lengths. Since £ c o tends to d s , while £ a bo diverges as £abo — £,abo/ (x — x u ) [20,31,33] we 
obtain 70 = d s /^ a bo, which is the critical amplitude of the anisotropy at the quantum superconductor to insulator 
transition at x u . The doping dependence of the relative isotope shifts is then traced back to a change of 70 and the 
shift Ax u of the underdoped limit and with that to a change of the mobile carrier concentration 8 = x — x u according 
to 
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A 7 (T e ) = A70 + Axu = A70 + Axu ^ 

T( Tc ) ^ 6 ^ x m (l±y/l- T c /T c (x m )) ' 

where we invoked the empirical relation between the hole concentration x and T c due to Presland et al. [34]. x m ~ 0.16 
denotes optimum doping. This leads to the important conclusion that the doping dependence of the isotope effects 
in the cuprates stem from the shift of the underdoped limit. Finally, combined with Eq.(5) we obtain 

A 7 (T C ) A 7o Ax u AT C A£ ab0 A£ c0 Af c0 , 2AA afc0 

— T^FV ~ 1 r - = 1 ~ a = — Tf~ - —<■ 7- — 7 1 T > \P) 

7l J cJ 70 S J c ?abO t,cO t,cO ^abO 

relating the various relative shifts to the scaling factor a. With our estimates AT C /T C ~ —0.014 (x — 0.15) and 
AT C /T C ~ —0.07 (x = 0.08) for La2- x Sr x Cu04 and relations (7) and (8) we obtain with x u — 0.05 for the essential, 
but material dependent parameters, determining the doping and T c dependence of ar c = — (M/ AM) AT C /T C the 
values A70/70 — —0.01 and Ax u ~ 0.0024. To illustrate this feature and to check the generic significance of this 
scenario further we consider the copper isotope effect on T c in YF^CuaOy-a. As aforementioned, our scaling analysis 
of the susceptibility data of Zhao et al. [4] reveals full consistency with 65 x (T) = 63 x ( a T) for all doping concentrations 
within experimental error. Since \ — m/H the implications are equivalent to those derived for the magnetization 
and Eq.(8) applies as well. The resulting estimates for ar c = — (M/AM)AT C /T C are included in Fig. 3 and compared 
with those obtained from the traditional extrapolation approach [4]. More importantly, given A70/70 and Ax u the T c 
dependence of ar c is readily calculated with the aid of Eqs.(7) and (8). As shown in Fig. 3 in terms of the dashed line, 
agreement is achieved with A70/70 — —0.0082 and Ax u ~ 0.0012. In comparison with Yi- y Pr y Ba2Cu3 16 ' 18 0y_5 
the data of Franck et al. [35] yields A70/70 ^ -0.0060 and Ax u ~ 0.0019. 
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FIG. 3. 7 (T c ) vs. x for I^-^Sr^CuCU taken from [18,26-30] (•) and a Tc vs. T c /T c (x m ) for YBa 2 63 ' 65 Cu 3 7 -5 taken from 
Zhao et al. [4](B). The dash-dot line is Eq.(6) with x u = 0.05 and 70 = 2. The arrow indicates the flow to the superconductor 
to insulator transition. The open squares result from the scaling analysis of the susceptibility data for the samples with 
7-5 = 6.94, 6.75, 6.63 and 6.48 in terms of 65 x(T) = 63 x (aT) yielding a ~ 1.0006, 0.994, 0.990 and 0.988, respectively. The 
dashed and dotted curves result from Eqs.(7) and (8) with A70/70 = -0.0082, Ax u = 0.0012 and T c (x m ) = 92.37 K. 

In summary we have seen that near T c , where 3D-XY fluctuations are essential, the isotope effects on various 
critical properties are not independent but related by universal relations. Together with the observation, that data 
taken at fixed magnetic field and on samples with 16 or 18 0, as well as with 63 Cu or 65 Cu, collapse near criticality 
within experimental error, when the temperature is rescaled appropriately, we derived an additional relationship. It 
reveals the essential relevance of the anisotropy 7. Indeed, the relative shift of T c was found to mirror that of the 
anisotropy 7. As a consequence, the temperature shift of the magnetization, susceptibility and the magnetic torque at 
fixed magnetic field does not provide estimates for the change of the transition temperature only, as hitherto assumed 
[1-9]. Together with the generic behavior of the anisotropy [18,26-30], the doping dependence of the isotope effects 
was traced back to a change of the underdoped limit Ax u , or in other words, to a shift of the phase diagram in the 
temperature-doping plane towards a slightly higher dopant concentration, along with a reduction of the mobile charge 
carrier concentration. This contribution leads to a negative shift of T c . We identified a positive shift as well. It stems 
from the change of the critical amplitude of the anisotropy 70 at the quantum superconductor to insulator transition. 
The magnitude and proportion of these contributions is controlled by Ax u and A70. Their values turned out to 
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be material dependent. In any case, they control the isotope effects and remain to be understood microscopically. 
However, the emerging essential role of the anisotropy represents a serious problem for two dimensional models as 
candidates to explain superconductivity in the cuprates, and serves as a constraint on future work towards a complete 
understanding. In addition, isotope exchange leads unavoidably to lattice distortions. Their coupling with the in-plane 
penetration depth was recently established [24]. 
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